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Abstract. We present a study of the entanglement properties of Gaussian cluster 
states, proposed as a universal resource for continuous- variable quantum computing. 
A central aim is to compare mathematically-idealized cluster states defined using 
quadrature eigenstates, which have infinite squeezing and cannot exist in nature, with 
Gaussian approximations which are experimentally accessible. Adopting widely-used 
definitions, we first review the key concepts, by analysing a process of teleportation 
along a continuous-variable quantum wire in the language of matrix product states. 
Next we consider the bipartite entanglement properties of the wire, providing analytic 
results. We proceed to grid cluster states, which are universal for the qubit case. To 
extend our analysis of the bipartite entanglement, we adopt the entropic-entanglement 
width, a specialized entanglement measure introduced recently by Van den Nest M 
et al. 2006 Phys. Rev. Lett. 97 150504, adapting their definition to the continuous- 
variable context. Finally we add the effects of photonic loss, extending our arguments 
to mixed states. Cumulatively our results point to key differences in the properties of 
idealized and Gaussian cluster states. Even modest loss rates are found to strongly 
limit the amount of entanglement. We discuss the implications for the potential of 
continuous-variable analogues of measurement-based quantum computation. 



PACS numbers: 03.67.Lx, 03.65.Ud, 42.50.-p, 42.50.Ex 
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1. Introduction 

Continuous-variable (CV) quantum computing offers a promising route to optical 
quantum computing, being naturally suited to common experimental techniques using 
squeezed light and homodyne detection [Tj. CV Gaussian cluster states have been 
proposed as a resource for universal quantum computation using the one-way paradigm 
[21 [3] . Several theoretical proposals address the problems of efficient generation of these 
cluster states as well as performing the necessary local operations [HE], and some proof- 
of-principle experiments have been reported [6] . These constitute part of a programme of 
adapting techniques from the context of discrete-variable (DV) - primarily qubit-based 
- quantum computing to phase space defined for infinite-dimensional modes, examples 
including quantum teleportation and other standard quantum-information protocols, 
and the matrix product state (MPS) representation [?]. 

Often concepts are generalized from the qubit case to phase space by reference to 
some limit of infinite squeezing. Although such a mathematical procedure is intuitively 
appealing, the arguments typically used are not exact for physical states, since finitely- 
squeezed approximations of non-physical infinitely-squeezed states may demonstrate 
behaviour departing rapidly from the ideal. For example, this can lead to cumulative 
errors in a computation (8j. Strong squeezing is difficult to achieve in practice, although 
a high-squeezing factor was reported recently in [9]. Squeezed states can be degraded 
by photonic loss and are commonly considered to be fragile, although this depends 
on the specific task and figure of merit [10J. Together these issues place demanding 
requirements on error mitigation and fault tolerance strategies for implementing CV 
quantum computation. An important line of investigation is reevaluate the usefulness of 
Gaussian cluster states as resource states when realistic assumptions on finite squeezing 
are applied. To this end, the scaling of localisable entanglement was considered in [8], 
the local complementation rule in pjj and a CV implementation of the Deutsch-Jozsa 
algorithm in [T2] . 

In this work, we consider several bipartite entanglement properties of Gaussian 
cluster states, making comparisons with their discrete analogues. As an application, 
we consider the recently-introduced entropic-entanglement width (EW) measure, which 
provides a necessary condition for universality for families of resource states p2]. We 
begin in section [2] by introducing idealized CV cluster states, and their approximations 
using Gaussian states. An intuitive demonstration of the problems of finite squeezing is 
presented by looking at a process of teleportation. In section [3j we look at the entropic 
entanglement for all possible bipartitions of a quantum wire, focusing on the case of 
pure states. The results of this are then applied, in section |1J to approximating the EW 
for a Gaussian grid cluster state, the qubit analogue of which is a universal resource 
for DV quantum computation. In section [5] we present numerical results on the effects 
of photon loss on the bipartite entanglement, adopting the logarithmic negativity as 
an entanglement monotone suitable for the mixed-state case. Finally, in section [6] we 
discuss the relevance of these results for the prospects of scalable CV measurement-based 



Bipartite Entanglement in Continuous-Variable Cluster States 



3 



quantum computation. 



2. Idealized CV and Gaussian cluster states and the effects of finite 
squeezing on coherent information transport 

Let us first recall a widely-used definition of (idealized) CV cluster states |3J 4]||] For 
each mode labelled i (termed here a "qumode"), with annihilation operator Sj, we work 
in a convention for which di,a\ = 1, Xi = (a\ + a,i j / \/2 and Pi = i (a\ — a^j / y/2. Given 
a graph with N vertices, for which the neighbourhood of vertex r is denoted n(r), the 
corresponding CV cluster state is, 

N 

|CVcluster) 00 = n J] \p = 0), ■ ■ • \p = 0) N , (1) 

r=l s£n(r)\s>r 

where the basis states are momentum-quadrature eigenstates satisfying P\p)=p\p) (for 
further explanation see [L7]), an d CV controlled-Z operations acting on pairs of modes 
are defined by Uq^ =exp (iX T X s j . 

To explain the analogy with conventional quantum information processing using 
qubit cluster states, we consider as an example a process of "one-bit teleportation" 
[18] along a computational quantum wire, and apply the MPS representation as in [19] . 
The MPS formalism provides, for any choice of physical computational device, a natural 
representation of the information-processing implemented by local measurements acting 
on the underlying computational resource states. Accordingly, a qubit wire of iV sites, 
with initial input state \ipm) (possibly unknown) defines a corresponding "correlation 
space" and matrices A^' 1 '^-] acting on it as: 

N 

h/W =^" 1) ---^Z 1) |^in) 1 (g)|+), 

i=2 

N 

= J] (Wl^^ N Wi) ® kA. (2) 

where the x r label the computational basis states with values and 1, |±) = 
(|0) ± |1)) /y/2, and controlled-Z operations are defined by U C z = |00)(00| + |01)(01| + 
|10)(10|-|11)(11|. Writing U$£ 1,r) \+} r+1 = |0) p A^ r+1 ^ [0] + |l) r A^ r+1 ^ [1], we see that 
MPS matrices acting between sites r and r+1 are given by A^ r+1 ' r ^ [0] = |+) r+1 (0| r and 

A^>r)[l] = |-} r+1 (l[ r . 

Suppose we choose another basis labelled {|</v)} and perform a measurement in 
this basis on site r, leading to projection \<f) r ){4>r\r, the corresponding MPS matrix is 
given by A^ r+1 ' r ^ [<f> r ] = ^2 X {4>r\ r \ x r) r A < - r+1 ' r ) [x r \ (as follows directly from the definition 
of correlation space given above.) To implement the paradigm of measurement-based 

| An alternative approach might use a qubit encoding of computational states into the infinite- 
dimensional modes, as discussed in |15j . An implementation using common-optical techniques was 
recently proposed (TB]. We do not pursue such encoding schemes here. 
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quantum computing, we need to choose a measurement basis so that the A^ -1-1 '^^] are 
unitary (for example A[0, 1] are only rank 1). Choosing the customary measurement 
basis {|±)} we find: 

^ (r+1 ' r) [+,] = (l+U(o| r + l-> r+1 (i| r )/v^ 

A (r+1 > r) [-,] = (|+> r+1 (0| r - |-> r+1 (l| r ) /V2. (3) 

These matrices are unitary (up to a normalization factor), involving a change-of-basis 
transformation (a Hadamard matrix) and, in the second case, an additional initial Pauli- 
Z operation. Hence measurements in the {|±)| basis drive unitary transport of 
along the wire, and the randomness of the measurement outcomes can be corrected by 
a local unitary operation at every step (thus implementing deterministic computation 
using feedforward control). Measurements along other axes in the x-y plane of the Bloch 
sphere allow one to generate any SU(2) logical gate, from which a fully universal set can 
be achieved by addition of an entangling gate. 

We proceed now to make the analogy for transporting an initial unknown "qumode" 
state \ip in ) along an idealized CV quantum wire, using the language of generalized 
functions. Defining the wire using ([I]), we identify the computational basis states \x r ) r 
as position-quadrature eigenstates (for which X\X = x) — x\X — x), the x-labels are now 
continuous, (X = x\X = y) = 5 (x— y) and dx \X—x)(X — x\ — I). The momentum 
quadrature (generalized) basis now plays the role of the |±) basis for the qubit case, 
and the quadrature bases are related by (X = x\P = p) = exp (ixp) /y/2n. Applying 
the MPS formalism for the CV case, the A[-} matrices in the computational basis are 
given by: A^ r+1,r ^[x r ] = |P = x r ), r+1 (X = x r \ r . For measurements in the momentum basis 
(corresponding to homodyne detection): 

At +1 ' r) i-Pr] = (2tt)- 1/2 I J dx r \P = x r ) r+1 (X = x r \\ exp (ip r X^ , (4) 

which is unitary up to a normalization factor. Analogy is achieved in the ideal case with 
the matrices A[±] for the qubit case: exp (^isX^j induces a translation \P) \— > \P+s) and 

plays a similar role to Pauli- Z (likewise exp(-itP) induces a translation \X) i-> \X+t)) 
and can be identified with Pauli-A); the Fourier transform operation corresponds to the 
"always-on" Hadamard operator. 

Going further, the Weyl-Heisenberg group of displacement operators (acting on 
a qumode) might be identified with the Pauli group (acting on a qubit). Gaussian 
operations, acting by conjugacy, transform each quadrature operator to a linear 
combination of the quadrature operators (possibly with an additional translation and 
all coefficients being real). They therefore transform every displacement operator to 
another displacement operator, and stand in the same relation to the Weyl-Heisenberg 
group as the Clifford group does to the Pauli group [20J. A non-Gaussian interaction 
characterized by a third-order Hamiltonian, or a process of photocounting jl], can 
achieve universal CV quantum computation in the sense set out in pQ, in a similar 
way as a non-Clifford operation such as a single-qubit 7r/8-phase gate is required for 
universal quantum computing with qubit systems. 
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Having made these analogies, a more realistic assumption is to replace the 
quadrature eigenstates for the CV cluster state with Gaussian states having finite 
squeezing. (We retain projections in the momentum-quadrature basis for the detection 
component, as these can be realised to good approximation by homodyne detection.) 
For each mode in ([!]) replace |P = 0) h-> (£) |vac), where the unitary one-mode 
squeezing operator is defined by (£) = exp [(— £a^ 2 + (*a 2 ) /2] and C— \C\exp(i(p). 
The following relations encapsulate the effect of S^(() in phase space. 

S^XS {1) = [cosh (ICI) -cos (<p) sinh (|C|)] X - sin (<p) sinh (|C|) P 
SWipsW = - S in (ip) sinh (|C|) X + [cosh (|C|) + cos (<p) sinh (|C|)] P. (5) 

Setting C = — | C| achieves a squeezing for momentum, and an anti-squeezing for 
position. The squeezed vacuum has only even photon number components, and hence 
the expectation value for the quadrature operator in any phase-space direction is zero. 

It is now possible to write down directly MPS matrices A[-} corresponding to 
Q above, for physical Gaussian cluster states. Expanding the squeezed vacuum 
in terms of momentum-quadrature eigenstates, dp\P = p)(P—p\ | — |C|) = 



TOO 
-DC 

[27T (A 2 p)] dp exp — 4( ^ 2 s \P = p), having a Gaussian distribution about P = 



with A 2 p = [exp (-2 |C|)] /2. Then, 

1 physical 



^SbiihPr] cx <! / dp exp 



exp (ipXr+l)\ A&W[-Pr].(6) 



4(A 2 p) 

We see that the effect of finite squeezing is to introduce an additional smearing in the 
computational basis (which has been mapped to the momentum basis after the action 
of the measurement on qumode r). The A[-] matrix no longer acts as a unitary map 
along the wire, and the errors arising from measurement are not simply correctable by 
feedforward. Every attempt to transport the original qumode another step along the 
Gaussian wire will introduce a further spreading. Hence the effect of finite squeezing is 
very detrimental to accurate coherent information transport, a point previously made in 
[H IHJ. We now proceed to investigate the computation power of Gaussian cluster states 
from another angle - namely their bipartite entanglement properties. 



3. Pure-state entanglement for bipartitions of the quantum wire 

In this section, we begin our study of the bipartite entanglement properties of 
Gaussian cluster states, disregarding decoherence processes for the time being. Since 
these quantum states are pure, we can adopt the entropic entanglement (EE) as 
a convenient measure; it is defined by EE(/3i,/32) — — tr [p re d log (p re d)] [|J where 
0i and /3 2 denote a bipartitioning of the vertices, and p ro d denotes either of the 
corresponding reduced density matrices. Entanglement monotones are invariant under 
local unitary transformations, and so controlled-Z operations within fli or 02 can safely 



§ Following common conventions for EE, we use the natural logarithm for CV calculations, and base 
2 for qubit calculations. 
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be disregarded, along with qumode/qubit vertices which become disconnected this 
way. Our basis for comparison is the EE properties for qubit cluster states, which 
can be determined for any given bipartition by simple counting arguments [H]. More 
specifically, let A denote the adjacency matrix of some finite graph (such that Aij is 1 
if vertices i and j are joined by an edge and otherwise), and A' denote the submatrix 
obtained by deleting from A rows corresponding to the elements of Pi and columns 
corresponding to /3 2 . A' then encodes information about edges between (3i and /3 2 , and 
EE((3i, P2) is given by the binary matrix rank of A' (i.e. arithmetic modulo 2, and where 
base-2 is used for the definition of EE) |21j . In a qubit cluster state, the contribution 
to EE{fix-i fi-i) from a particular vertex cannot be more than 1 whether it is singly or 
multiply bonded, since the maximum entropy of a qubit is 1. In the CV case, however, 
there is no upper bound to the entropy of an individual qumode and the entropy of 
entanglement can exhibit a richer range of behaviour, as we shall see below. 

For our CV calculations, we use phase-space methods [22J, defined for N modes as 
follows. The 2N canonical-coordinate operators and variables are denoted by column 
vectors 6 = (Xr, X 2 ; ■ ■ ■ A; A • • •) and £ = (X x ; X 2 ; ■ ■ • Pf, Pi • • •)■ Matrix £ is defined 

by Oj,Ok = iT,jk, where X, P = i. A (pure or mixed) state p is called Gaussian if 

its characteristic function, Xp (0 = t r pexp ^— i^'Soj , is completely determined by 
its first and second moments and can be written in terms of a displacement vector d 
and covariance matrix T, xp (0 =exp (— |£*£r£*£ + id*££), where tr [pO^j = d a , and 



ab 



2tr 



p (6 a - d a j (6 b - d t 



— iT, ab . A further condition applies: T represents a 

physical state if and only if T ± iT, > 0. Given a unitary transformation U on the state 
space for which U^OiU = Y^j SijOj, the matrix S is a real and symplectic (namely it 
satisfies the condition S t T l S = X). The transformation of a Gaussian state p 1— > UpU^ 
is equivalent to the mappings T i-> SYS* and d (->• Sd. 

For a Gaussian iV-mode state, the EE can be evaluated directly from the covariance 
matrix. Given a bipartition of the modes, the reduced state for one component is 
Gaussian, and the corresponding reduced covariance matrix r re( j is given by removing 
rows and columns corresponding to the modes which are traced out. Applying 
Williamson's theorem, which provides a normal-mode decomposition for real, symmetric 
and positive-definite matrices [23J, any covariance matrix T can be diagonalised via a 
symplectic transformation. The corresponding "symplectic eigenvalues" come in same- 
value pairs with positive values Ai, •••,Ajv (where Aj > 1 for physical states). The 
matrix Er has eigenvalues ±zAi,- • -,±iAjv, and can be used to compute the symplectic 
spectrum. The diagonalising symplectic transformation in phase-space corresponds to a 
unitary transformation on the state-space, relating the Gaussian state to thermal states 
in iV independent modes. These results allow a simple formula for EE to be derived, 
from the symplectic spectrum A£. ed of r re( j: 



EE = Y J 



k=l 



(7) 
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We now focus on the bipartite entanglement behaviour of a quantum wire. For any 
given bipartition, we disregard controlled-Z operations which do not contribute to the 
EE, and thereby break the wire into disconnected lengths of "zigzig-type" strings, for 
which alternate vertices switch bipartition component. The simplest zigzag has only 
two vertices and a single bond. Disconnected zigzags contribute additively to the total 
EE. Hence the EE for all possible bipartitions can be understood by considering only 
zigzag bipartitions. For a continuous zigzag bipartition of a iV qubit wire, the EE is 
easily seen to have value (N — l)/2 when N is odd, and N/2 when N is even. We 
progress now to the case of Gaussian quantum wires defined as section |2j We note that 
a related investigation is reported in [21] , which looks at the ground and thermal states 
of a closed harmonic chain with nearest-neighbour couplings. 

We first transform the iV-qumode Gaussian quantum wire by local unitary 
squeezing processes, absorbing the initial momentum-squeezing of each qumode into 
the CV controlled-Z operations: 

Sl n |G - wire) = J] (s^^aii) |0> , (8) 

where S a u = ®^f =1 S (— \(\) r and S^Ug® = exp (ie'^XaX^j . In other words, we 
can replace the single-mode squeezing and controlled-Z operators by a single two-mode 
squeezing operator, whose strength is given by the squeezing parameter B = e~~ 2 \ 
solely dependent upon the degree of squeezing of the initial single-mode states. By 
the Baker-Hausdorff lemma, the symplectic matrix for two arbitrary modes labelled i 
and j corresponding to gate operation U^'^(B) 



S(B) 



( 


1 








o\ 







1 













B 


1 





\ 


B 











(9) 



Let us begin by studying the simplest zigzag, namely two modes connected (across 
the bipartition) by a single bond. For state exp (iBXiX-zj 1 0) 1 1 0) 2 , the symplectic 

eigenvalue for the reduced density matrices is A = vT+ B 2 , and there is one contribution 
to the EE in (J7|). Provided A is not close to 1, a convenient approximate form for the 
entanglement entropy of this two-mode state (which upper bounds the exact expression, 
and converges rapidly to it with increasing A) can be readily derived, 

~ In f £ J + 1. (10) 



2, 

Hence, in terms of the squeezing parameter B, the entanglement entropy for B ^> 1 has 
the simple form EE S ~ \n(B/2) + 1. 
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We can calculate the entanglement entropy for certain related graph states, 
across appropriate bipartitions, by reducing them to the single bond case, but with 
a modified squeezing parameter. Consider a star Gaussian cluster state with N — 1 
edges: exp (iBX x X^ exp (iBXxX^ • ■ ■ exp (iBXtX^ |0)i|0) 2 ■ ■ ■ |0>jv (1 labels the 
centre mode). Taking the bipartition of mode 1 versus all other modes, as illustrated 
in figure f]Ta), we can transform this state via a local unitary (with respect to the 
bipartition) to a two-mode single bond state. We achieve this via the mapping 
(X 2 + X 3 ---Xn) i — y \/N — \X' 2 . This can be achieved using a mode unitary 

beam-splitter network implementing the finite Fourier transform (defined by the linear 
transformation a s h-> J2t F sta t with (F) st = [exp(i27T (s- 1) (t-1) /(N-1))] /y/N-1) [25]. 
Hence we see that the star state has the same EE as the single bond state with squeezing 
enhanced by a factor of \/N — 1. For initial squeezing B, the EE for the star state is 
approximately EE stai ~ In (y/N-lB/2) + 1. 

Another example is the completely-connected graph state, defined by 

n^j ex P (iBXiXj) 1 0) 1 1 0} 2 ■ ■ ■ |0)jv, where a bonding operation has been applied 
once to every pair of modes. Given a bipartition of the modes into components with 
size r and N—r, as illustrated in figure[ljb), we can use a similar trick as above to map the 
state into a single bonded pair with increased bond strength B \— > ^/r\jN — rB. The en- 
tanglement entropy in the strongly-squeezed limit is EE comp \ ete ~ \ia(y/ry/N —rB/2) + 1. 
A maximizing bipartition then has components of equal size. We have found that the 
scaling for the EE is logarithmic in the system size, both for the star and completely- 
connected configurations of the Gaussian cluster state. Following similar arguments to 
|14j . we can immediately conclude that these two families of states do not constitute 
universal resources, for any value of the squeezing parameter B. 

The qubit analogies for these states behave rather differently. In particular, the 
iV-qubit cluster states in the star and completely-connected configurations share the 
same entanglement properties. This is true since a star graph can be converted to the 
completely-connected graph by the local complementation rule (applied to the centre 
qubit), which implies the existence of a local Clifford operation transforming one into 
the other J2JJ. Therefore, for every nontrivial bipartition of these qubit states, the EE 
has value 1. 

We will now consider a wire of length N, initially in the vacuum state, subject 
to operations exp yiBXjXj+xj on pairs of qumodes (1, 2), (2, 3),- • -,(N — l,N). We 
bipartition the wire as a zigzag with all odd-numbered modes constituting one half of 
the chosen bipartition (as illustrated in figure []Jc)). The reduced covariance matrix for 
the even-numbered modes is of the block form, 




where the blocks have dimension iVo = by when iV is odd, and = | by | 
when N is even. denotes the identity matrix, and T/v is a tridiagonal matrix. When 
iV is odd, the leading diagonal of T No has constant value (1 + 2_B 2 ), and the adjacent 
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Figure 1. Three classes of iV-qumode Gaussian cluster states, generated by CV 
controllcd-Z operations with parameter B. (a) Star configuration with bipartition 
defined by central qumode 1; (b) completely-connected configuration bipartitioned 
with components of size r and N—r (controllcd-Z operators that do not contribute to 
the EE have been omitted); (c) zigzag configuration bipartitioned into qumodcs with 
odd or even labels. For calculating the EE, cases (a) and (b) are equivalent to the 
bipartite entanglement for a pair of qumodes acted upon by a controlled-Z operation 
with B H> ^N — IB and B h-> ^r^fN—rB respectively. 



diagonals have value B 2 everywhere. The N even case is similar, with the bottom right 
corner element altered to (1 + B 2 ). For the symplectic spectrum corresponding to the 
zigzag bipartition, we need the eigenvalues A k of Sr^g zag (N). These also solve the 
eigenvalue equation Tn q v = —K\v. The spectrum of T/v is easily computed for odd and 
even N cases |2"E1. and we find, 



At 




2B 2 1 



cos 



2k 



7T 



N+l 



(12) 



where k = 1 • • ■ No; the EE can be computed from the absolute values using ([7]). This 
expression appears rather complicated, but we can derive a much simpler form valid in 
the limit of large squeezing and many modes. We shall focus only on the case where 
TV is odd since it will play an important role in calculating the entropic-entanglement 
width for a square-lattice Gaussian cluster state later in this article. 

We begin by taking the high-squeezing approximation (10) for every term in ([7]). 



iV 



EE, 



zigzag 

k=l 

45> 



k=l 



1 + 2B Z 1+cos 



2kn 
N+l 



+ jV [l-ln(2)]. 

In the limit that Nq is large, we can approximate this sum by an integral to obtain, 



EE, 



zigzag 



(N,B) 



No 
2vr 



ln[l + 2B 2 (1 + cos(x))]dx + jV [l-ln(2)] 
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Figure 2. Two examples of decompositions, each denned by a subcubic tree (for 
which every vertex has one or three incident edges) and a bijection from the tree 
leaves to the qumodes. Removing an edge defines a bipartition. Vertices with only 
two incident edges are associated with only one bipartition, and can be removed by 
merging the edges. Although not necessary for the definition of the EW, subcubic 
trees can be illustrated as rooted trees, for which an arbitrary non-leaf vertex has been 
designated the root, defining a partial ordering. 



and upon integrating [27], 

EE zigzag (N,B)~N 



and hence, 



EE zigzag (N,B) « N 



l--ln(2) + ln ( \/l + 2B 2 + Vl+4B 2 



In \ 1 + 2B 2 + Vl+4:B 2 



(13) 



(14) 



We find that (14) converges quickly to the value of the exact expression as B and N 
grow large. 



4. Entropic-entanglement width for the grid state 

In this section we consider the entanglement properties of the family of Gaussian cluster 
states in the Z-by-Z-grid configuration, which is well known to be universal for DV 
quantum computing. To do this we adopt the entropic-entanglement width (EW), 
introduced by Van den Nest et al. [13] . The previous authors focused on the implications 
of the scaling of the EW, and provided a no-go result for (efficient) universality when 
the EW is bounded (or scales at most logarithmically) for a family of resource states. In 
this work we extend the definition of EW, replacing qubits with qumodes, and we look 
primarily at qualitative differences in the behaviour of Gaussian cluster states compared 
to qubit cluster states. 

To begin, we introduce the idea of using a subcubic tree (a tree graph having one 
or three edges incident at each vertex) to induce a set of bipartitions of a collection of 
qumode (or qubit) labelsjjj] The situation is illustrated in figure[2] A tree is an undirected 
graph in which any two vertices are connected by exactly one simple path. It follows 
from this definition that any tree is connected, without cycles, and that the number of 
edges is one less than the number of vertices. Vertices with only one incident edge are 
termed leaves. Each leaf is identified with one qumode label, and a tree together with 



The use of trees here is entirely unrelated to any graph appearing in the definition of the state. 
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Figure 3. (i) A construction for building a decomposition from a chosen partition of 
graph vertices. A simple tree is constructed for each partition subset, and leaves are 
identified with the elements. These trees are then "strung together" (suitably ordered) 
to form a subcubic tree and a decomposition for the entire graph. For qubit cluster 
states defined on the l-by-l grid, the width for the decomposition defined by diagonal 
subsets (ii) is l — l, and the decomposition is optimal. For the rectangular case (iii), 
the width for the decomposition is I. For Gaussian cluster states, the widths for (ii) 
and (iii) depend on the squeezing parameter. 



this labelling is termed a decomposition. Given a decomposition T for a iV-qumode 
state \ip), removal of an edge e disconnects the tree, and bipartitions the set of qumodes 
Q into subsets (3 e and Q \ (3 e . Each edge of T therefore defines a value for the EE 
across the corresponding bipartition. The maximum of these values for T is termed 
the width of T. The EW of \ip) is then defined as the minimum width for all possible 
decompositions: 

EW{\f)) = mm (waxEE(P e ,Q\ • (15) 

We thus see that any valid decomposition can provide an upper bound to the 
EW, but finding the optimal decomposition (and hence a tight value for the EW) 
can be difficult. For qubit cluster states, the EE has the special form explained in 
section |3j and hence the EW [H] coincides with a graph width parameter called the 
rankwidth [28]. For a one-dimensional cluster state (a "wire") the rankwidth is easily 
seen to have value 1. For the /-by-/ grid, by consideration of the bipartition defined 
by the main diagonal, one can show that an upper bound to the rankwidth is l — l. 
Recently, by a series of clever arguments which exploit the symmetries of the square 
grid structure, it has been shown that this bound is tight and that the rankwidth of 
the grid is / — 1 [29] . To evaluate the difficulty of computing exhaustively the EW for 
an arbitrary iV-qumode state, we count the total number of possible decompositions 
(for N = 2,3 there is only one possibility). Suppose that T^-i is a decomposition 
for (N — l) qumodes; it can be turned into a decomposition for iV qumodes by the 
following step: split an edge and insert a non-leaf vertex together with an adjacent leaf 
labelled N. This step adds two vertices and two edges. Furthermore, every edge of 
every (N— l)-qumode decomposition leads to a different iV-qumode decomposition, and 
all possible iV-qumode decompositions are generated by this method. It follows that 
there are (2JV-5) (2JV-7) • • • 1 = (2JV-5)!/ [2^ (iV-3)!] total decompositions for a 
iV-qumode state (N > 3), and direct evaluation of the EW is unfeasible unless N is 
small. 
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To address the problem of computing the EW for large grids, we select two candidate 
decompositions whose properties, in the case of qubit states, make them good candidates 
for optimal decompositions. Calculating the width of these decompositions provides 
upper bounds for the EW for the states. For small grids, we can then check the 
optimality of the decomposition via numerics. 

The construction we use is illustrated in figure [3](i)|l| Given a graph, we first 
partition the set of vertices Q into subsets (3i,---,(3 p (possibly of different sizes). 
For each = {vn, v^, • • •}, we associate a simple tree. Removing edges for this 
tree generates bipartitions ({vn},Q \ {%l})> ({ v n, v i2},Q \ {^1,^2}), and so on. We 
can then "string" these simple trees together in order, to yield a sub cubic tree for 
the whole of Q. Removing edges from the string generates bipartitions (Bi,Q \ 
Bi), (Bx U B 2 , Q \ (Bx U B 2 )), and so on. The two decompositions we choose use 
diagonal bipartitions (figure [3j^ii) ) and rectangular bipartitions (figure [3jiii) ) , where the 
bipartitions are strung together from top left to bottom right, and from top to bottom 
respectively. 

For the diagonal-type decomposition, the bipartition defined by the central diagonal 
is maximizing and determines the width of the decomposition (although it is not in 
general the only maximizing partition). In both DV and CV cases, by applying local 
unitaries within the bipartitions, the calculation of the EE reduces to that for the zigzag- 
type bipartition, whose properties for Gaussian cluster states we analysed in section [3j 
For the case of a qubit cluster state, the width of the diagonal-type decomposition 
takes the optimal value, I — 1. For the rectangle-type decomposition, the maximizing 
bipartition will be given by a row of vertices away from the boundary of the grid. 
Applying unitaries local with respect to the bi-partition we can remove all entangling 
"bonds" leaving / separate and unentangled 3-mode wire states. For the case of a qubit 
cluster state, the width of the rectangle-type decomposition is immediately seen to be 
/, and the decomposition is suboptimal. 

To calculate the width of these decompositions for a Gaussian cluster state of 
bond-strength B defined on the /-by-/ grid, we can directly use results derived in 
section [3j First we shall consider the diagonal-type decomposition. Via application of 
local unitaries with respect to the bipartition defined by the central diagonal, the state 
can be reduced to a (/ — 1) length zigzag cluster state with bond strength B' = y/2B. 



Then, assuming B is not small, we can directly use the approximate expressions (13) or 



more compactly (14) to achieve 



EE diag (l, B) « (l-l) In ^1+452 + ^1+85^ . (16) 

For the rectangle-type decomposition, by the application of further local unitaries 
(as for the star graph in the section [3]) we can reduce / length-3 zigzags to / singly- 
bonded pairs across the bipartition with bond strength y/2B. The corresponding EE is 



therefore / times that of a simple pair with this bond strength. We use (10) with the 
% Vit Jch'nck, private communication. 
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single symplectic eigenvalue A = y/l+2B 2 to directly obtain the following expression, a 
good approximation when A is not close to 1, 

EE Iect (l,B) ~ / [ln( v / l+2 J B 2 ) + l-ln(2) 

wZ In (Vl+2£ 2 ) . (17) 

For any given size of cluster state, we do not know whether either of these 
two decompositions is optimal, so we have investigated the behaviour for small grids 
numerically, evaluating the EW for the 3-by-3 and 4-by-4 grids. To do this, we use an 
equivalent definition of the EW based on recursion [30J. Letting Q denote the set of 
qumode labels, define the function ty(-) on all subsets X of Q as: 

w(X) = mm Yc x\Y^9,Y^x 

m&x{EE(Y, Q \ Y), EE{X \Y,Q\(X\Y)), w(Y),w(X \ Y)} 
when|X| > 2, 

EE(X, Q\X) when|X| = 1. (18) 

Similarly to the previous discussion on decompositions defined by subcubic trees, w(X) 
can be understood as the minimum width of a binary tree on the subset X of Q. A binary 
tree has one or three incident edges at every vertex other than the root; removing an edge 
from the binary tree defines a subset given by the leaves which are the descents of the 
edge. The recursive definition here works by combining binary trees which are already 
optimal on smaller subsets. The starting point is the values of the EE corresponding to 
singleton bipartitions of Q. For the case of X — Q, w(Q) now coincides with definition 
of the EW (an optimal binary tree here be converted into an optimal subcubic tree by 
removing the root and merging the two incident edges). ( 18 ) can be readily implemented 
using dynamic programming techniques [31] . The results are shown in figure |1J In both 
cases, we see that in the limit of high squeezing the EW tends to the diagonal-type 
decomposition. This behaviour is not surprising, since in the infinite-squeezing limit, the 
entanglement properties are qubit-like, and the diagonal decomposition is the optimum. 
We therefore conjecture that, in the limit of high squeezing, the EW of the a general 
grid will tend to that of the diagonal decomposition, namely (16). 

Given that very-high squeezing is currently not experimentally achievable, a more 
relevant limit is to keep squeezing held constant, and study the behaviour as grid 
size I gets very large. In this limit, we see a very different behaviour. For fixed 
B, there is always a grid-size V such that for all grid sizes I > V the width of the 
rectangular decomposition is smaller than the diagonal decomposition, and the diagonal 
decomposition can no longer be optimal. This behaviour can be seen in figure |5j The 
reason for this transition is that both EE^ iag (l,B) and EE Tect (l,B) are linear in I, but 
dEE^ ag (l, B)/dl > dEE Tect (l, B)/dl. Thus, for a given grid size, something akin to 
a phase transition occurs in the entanglement properties. We must emphasize though 
that we cannot know if the EE Tect (l, B) provides a tight upper bound to the EW in this 
limit. The true EW might be much lower. What this does indicate is that the EW for 
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Figure 4. Plots show values of the EW (computed numerically) for a range of 
the squeezing parameter B — exp (— 2£), compared to the width for the diagonal- 
type decomposition (blue) and the rectangular-type decomposition (red). The limit 
B — too, (C — > — oo), corresponds to perfect squeezing for the momentum quadratures 
for every qumode at the start. The limit B — >Q corresponds to perfect anti-squeezing 
in the momentum quadratures, and there is no entanglement in this limit. For 
both the 3-by-3 lattice, case (a), and the 4-by-4 lattice, case (b), the diagonal-type 
decomposition is optimal for larger values of B (as would be expected from the qubit 
case). However, for this to occur, the required squeezing must be greater for the larger 
lattice. 



squeezing B 




grid side I 



Figure 5. A comparison of width of the diagonal- type (blue) and rectangle- type (red) 
decompositions of Gaussian cluster states defined on Z-by-Z grids. These are determined 
respectively by the EE of a bipartition defined by a main diagonal, and by a full row 
(or column) away from the boundary. Each edge between qumodes, labelled i and j, 
corresponds to the CV controlled-Z operation exp(iBXiXj). As / increases with B 
fixed, the EE for the bipartition defined by the main diagonal becomes greater than 
that for a bipartition defined by a full row. In contrast, for the qubit case, the EE is 
always one less for the diagonal case. 



large, fixed-squeezing grids has a rather different behaviour to either the qubit cluster 
state or the infinitely-squeezed case. 
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5. Bipartite entanglement in the presence of photonic loss 

In this section, we consider how our previous results are affected by photonic loss, a 
principle source of error in optical experiments, especially affecting squeezing operations 
and storage (necessary for adaptive measurement). To begin, we consider how a loss 
channel acts on an arbitrary iV-mode Gaussian state with covariance matrix and 
displacement vector dj. We adopt the standard loss model, for which a mode r couples 
to an ancilla mode e, initially the vacuum, via a beam-splitter-type interaction. The 
beam-splitter transformation acts on the mode operators as a r y ^dr + ^l—rje, with 
transmissivity parameter rj. A partial trace is taken for mode e at the end. This loss 
channel is Gaussian (that is to say it transforms Gaussian states to Gaussian states), 
and we write T-° ss and d l ° ss for the covariance matrix and displacement vector for the 
state at the end. If independent losses are applied to every mode with the same value 
for rj, the following transformation is readily derived for the phase-space formalism: 

r l°ss =r]T + ( 1 _ 7] } J ; 

d loss = ^jd. (19) 

/ here denotes the identity matrix. The displacement vector can be disregarded, since 
it corresponds to a local property and plays no role in calculating the entanglement. 

The loss channel will be assumed to act on each mode of a Gaussian cluster 
state after it has been generated from initial squeezing operations on each qumode 
with parameter ( = — |C|, and CV controlled-Z operations between pairs of qumodes. 
Squeezing operations do not commute with the loss channel so, unlike the previous 
entanglement calculations, we can not simply "undo" entangling operations within 
bipartitions to simplify our calculations. Furthermore, arguments used above for 
disregarding CV controlled-Z operations within a chosen bipartition, and for combining 
multiple operations incident at a particular qumode, can no longer be assumed to apply. 
To investigate the bipartite entanglement for lossy-Gaussian cluster states, we need an 
entanglement monotone defined for mixed states. We choose the logarithmic negativity 
(LN), defined for a state p by LN = In (| |p pt | |i), where || • ||i denotes the trace norm 
and is equal to the sum of the absolute values of the eigenvalues, and pt denotes the 
partial transposition operation on either party [32]. As previously, we use the natural 
logarithm for the CV case. In operational terms, the LN provides an upper bound on 
the distillable entanglement. 

We can use the LN to define an entanglement width for mixed states. Applying 
the central ideas of Van den Nest et al. [13] , our new entanglement width LNW must 
satisfy the following condition: if state p with N modes can be converted to state p' 
with N' < N modes deterministically, by local operations and classical communication 
(LOCC), then LNW(p) > LNW(p'). This condition defines a "type-II entanglement 
monotone" , in the terminology of Van den Nest et al. Differently from those authors, 
we consider CV states defined on qumodes rather than DV states defined on qubits, 
and we go beyond the pure state case, defining our type-II monotone on mixed states. 
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However we retain the two crucial elements of type-II definition: the interconversion 
by LOCC is assumed to be deterministic and there is no averaging of probabilistic 
outcomes; states are compared on subsystems of possible different size i.e. the end state 
can be defined on fewer qumodes. To verify that in fact our definition of LNW does 
in fact satisfy the type-II definition, it is sufficient to check two facts. First, we need 
to the special case that LNW(p) > LNW(p') when N' = N. This follows immediately 
from the generic min-max construction used to define the entanglement width, since 
LN is an entanglement monotone in the conventional sense (specifically that it is non- 
increasing on average under non-deterministic LOCC) [33]. Second, we need to check 
that the LNW is invariant under the addition of an uncorrelated additional party i.e. 
LNW(p) ®LNW(p® |0)(0|jv+i)- This is easily verified since, for an arbitrary bipartition, 
the value of the LN is unchanged by the addition of an extra party. 
We thus define the logarithmic-negativity width (LNW) as follows: 



LNW{\p)) = mm (maxLN(p e , Q\f3 e 



(20) 



where the width of a decomposition T is defined as the maximum value of the LN for 
bipartitions of the set of qumode labels Q, each one defined by removing an edge e of 
T. Again the corresponding width, LNW, is defined as the minimum width over all 
possible decompositions. 

We now compute the LNW for a lossy- Gaussian cluster state in a Z-by-Z-grid 
configuration using the phase-space formalism. As previously, we use B = exp(— 2() 
as the relevant squeezing parameter. Letting T l °^ d denote the covariance matrix for the 
state, we write: 



■ploss 
grid 



7] ( 

_edges{ij} 
h (1-7J)/. 




c(m) 



cdgcs{i,j} 



(21) 



S^z^ denotes a CV controlled-Z operation with qumodes with labels i and j, and 
has the value 1 on the main diagonal and at (i + l 2 ,j) and (j + l 2 ,i). In the phase- 
space formalism, the partial transposition operation acts to change the sign of the 
corresponding momenta. Given a bipartition of the qumode labels (/3i, /3 2 ), the operator 
after partial transposition (over f3\ say) is of the same form as a Gaussian density matrix 



with covariance matrix T 



loss,PT 



Prg^P, where P is a diagonal matrix having entries 



grid 

— 1 at positions {(q + l 2 , q + l 2 )\q G f3\}, and value 1 for the remaining diagonal entries 



Williamson's theorem can be applied to provide a global symplectic diagonalisation 
1 ;,> " 4 ••--'»• >i- ■<••*•;<•• >'-'•*'"• x1 ' 1 Modes with symplectic 



of r^'^ 1 in terms of symplectic eigenvalues A^ T 



\PT 
A [2 ■ 



eigenvalues Af T > 1 have value 1 for the trace-norm and do not contribute to the 
LN. Cases with < Af T < 1 contradict the uncertainty relation AXAP > 1/2, and 
contribute l/Af T to the trace norm [32]. Overall: 

i 2 

LN = J2 max [- ln ( A f T ) > 0] • ( 22 ) 



i=l 
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Figure 6. The LNW is plotted (dots) versus the squeezing parameter B for the 
(a) 3-by-3 and (b) 4-by-4 grid configurations. Each subplot corresponds to a different 
value for transmissivity parameter 77 (77 = 1 is the lossless case) . Blue and red lines are 
the LN for the diagonal and rectangular decompositions respectively. 
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(i) transmissivity r) = 1 : (ii) transmissivity t] = 0.9: (iii) transmissivity r] = 0.5: 




grid side I " grid M'cfe / 15 grid side I 15 



Figure 7. The LN is plotted for the widths of diagonal (blue) and rectangular (red) 
decompositions for /-by-/ grid configurations and for a range of the squeezing parameter 
B. Each subplot corresponds to a different value for transmissivity parameter 77 (77 = 1 
is the lossless case). 

In figure |6j Gaussian cluster states in the 3-by-3 and 4-by-4 grid configurations are 
compared. When loss is absent or not too large, (i.e. transmissivity 77 = 0.9 or 0.5), 
the diagonal decomposition is optimal, or close-to-optimal across the range of squeezing 
investigated. The convergence of the width of the diagonal decomposition to the LNW 
takes slightly longer for the larger grid but still occurs rapidly. However, there is a 
large drop entanglement as loss is added — for example when 10% loss is added the 
LNW falls to roughly half the value of the non-dissipative case across the range of 
squeezing considered. For the case of extreme loss, for which transmissivity rj = 0.25, 
the diagonal and rectangular decompositions are seen to be suboptimal for larger values 
of the squeezing parameter B, indicating a change in entanglement structure compared 
to the idealized case (of infinitely-squeezed CV cluster states or qubit cluster states). 
It is seen also that the entanglement grows slowly for larger values of B when loss is 
added, compared to the non-dissipative case. This suggests that errors arising from 
finite squeezing cannot simply be countered by achieving strong squeezing at the start; 
the infinite-squeezing limit cannot be assumed to behave as in the non-dissipative case. 

Finally, the diagonal and rectangular decompositions are compared in figure [7] for 
l-by-l grid configurations up to I = 15. For the non-dissipative case (i), rj = 1, similar 
behaviour is seen using the LN as a measure of entanglement as was previously observed 
for the EE (figure [5]). Specifically, for a fixed value of the squeezed parameter B but 
increasing grid size, the diagonal decomposition eventually has more entanglement than 
the rectangular decomposition. The phenomena becomes less marked as loss is added 
(cases (ii) transmissivity 77 = 0.9 and (iii) 77 = 0.5). In the lossy cases, the key features 
are the large drop in entanglement (as observed previously for the 3-by-3 and 4-by-4 
grids), and the linear growth of the LNW with increasing values I (with B fixed). 
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6. Discussion and conclusions 

In summary, we have investigated the bipartite entanglement properties of Gaussian 
cluster states under finite squeezing and photonic loss, focusing on quantum wires and 
grid lattices. In the case of finite squeezing but no loss, we have derived analytic 
upper bounds for the entanglement (quantified using the EW), which we conjecture to 
be tight in the limit of high squeezing. We have also shown that a different limiting 
behaviour occurs in the limit of large grid size, when the initial per-site squeezing ( is 
held fixed, compared to the DV case. When loss is added to the picture, we have found 
numerically a large decrease in the available entanglement (quantified using the LNW 
for mixed states). Furthermore, we have provided numerical evidence that the LNW for 
a grid cluster state grows only very slowly as the squeezing parameter _B = exp(— 2() is 
increased to large values whenever losses are present. In all cases however, analytic and 
numeric evidence point to a linear increase in the entanglement for /-by-/ grid states 
with respect to the grid-side / (B held fixed). 

It is worth considering the implications of these results for the prospects of scalable 
measurement-based quantum computation (MBQC) using Gaussian cluster states as a 
resource. In [13], Van den Nest et al. identified a set of criteria which every DV state must 
satisfy in order to be a "universal state preparator" (USP) under local measurements. 
This requirement, which is stronger than the weakest definition of a universality, (which 
would only specify a quantum processor providing classical statistical data at the 
output), is nevertheless satisfied by most known variants of MBQC [31] • The criteria are 
derived from type-II entanglement monotones, entanglement measures which are non- 
increasing under deterministic LOCC, where subsystems of possibly-different system 
size are compared. Every type-II entanglement monotone provides a no-go theorem 
which can exclude certain families of states as efficient USPs (for example see theorem 
9 in [13]). Variants of the entanglement width, defined in this paper in terms of the 
entropy of entanglement and the logarithmic negativity, constitute type-II monotones, 
and we can apply similar arguments about USPs to Gaussian cluster states]^] 

The Van den Nest et al. criteria are derived by reference to a known USP resource 
state, typically a grid (or rectangular) cluster state. The essence of the argument can be 
paraphrased as follows (for a formal description see [13] )■ An efficient USP must be able 
to generate (efficiently) any state creatable (efficiently) via a quantum circuit. Since this 
includes the family of grid cluster states itself, this implies necessary criteria relating 
to the scaling of their entanglement properties. For example, the EW for any efficient 
USP resource state must scale faster than logarithmically in the system size. We would 
wish to apply a similar criteria to Gaussian cluster states. A first choice which one must 
make is whether to compare these potential resource states with a known universal qubit 
resource (for example the same grid cluster states) or a suitable CV state. Choosing an 
appropriate CV state is problematic since, so far, the only generally-accepted USP for 
the CV case is the infinitely-squeezed cluster state. This state (and the high-squeezing 

+ See [8] for an analysis of the Gaussian localisable entanglement in these states. 
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limit of finitely-squeezed Gaussian states) cannot be used for comparison, since the 
entanglement can be unbounded even for the two-mode casej^] Instead, we shall make 
our comparison with qubit cluster states, and use directly the criteria derived in jT3] , 
namely that the entanglement width must be unbounded (for universality) and must 
scale faster than logarithmically (for efficient universality). 

For the pure-Gaussian cluster states the diagonal and rectangular decompositions 
both provide an upper bound to the entanglement width. As can be seen from (16) 
and (17) the entanglement in both these cases scales linearly in grid size, satisfying the 
criteria for (in principle) efficient conversion to qubit cluster states. These measures 
only provide upper bounds to the entanglement width, as we have not proven that 
these bounds are tight. They do, however, allow us to extract some indication of the 
capability of these states to support MBQC. In order to create a l-by-l qubit cluster 
state, the Gaussian resource state would have to possess an entanglement width of at 
least (I — 1) ln(2). This means that the prefactor in the linear EW scaling rules for the 
Gaussian cluster states can be thought of as a "best-case conversion rate" , quantifying 
the largest qubit cluster state which could be created via local measurements from a 
given Gaussian cluster state. For example, the central diagonal bipartition gives us, 
via (16), an upper bound on the conversion rate of In (y/l+4B 2 + >/l + 8B 2 } /ln(2). 
We find that a squeezing factor of B 0.54 (corresponding to some anti-squeezing of 
the momentum quadrature) provides a conversion rate of approximately unity. Note 
that even a state constructed from initially unsqueezed modes undergoing a full QND 
interaction would achieve this, and such states may be feasible to construct via a network 
of beam-splitters acting on single-mode squeezed states of currently achievable squeezing 
[35] . We emphasise, however, that this is a necessary condition only. It does not 
guarantee that such interconversion is possible, and certainly not that the measurements 
needed to do so could be achieved via linear optics. 

To conclude, while the results of our investigation have led us to less pessimistic 
conclusions than other recent investigations into Gaussian cluster state entanglement, 
these results taken together with [8] do suggest that utilising Gaussian cluster states for 
MBQC may require an approach radically different to the traditional qubit approaches 
[21 US]. Furthermore, they remind us that no scalable scheme for MBQC with finitely- 
squeezed cluster states has been proposed so far. We have seen for the entanglement 
width that the minimizing decomposition, in the limit of large grids, is not the same 
as the the optimal decomposition in the limit of infinite squeezing. This indicates 
that some aspects of the entanglement for finitely-squeezed states are qualitatively 
different to the infinitely-squeezed limit. The results of section [5] show us that the 
entanglement width degrades swiftly in the presence of loss, and that coding for tolerance 
of loss errors will need to be an important component of such schemes. On the other 
hand, our results show that, with only a modest initial squeezing, it may be possible 

* A better family of comparison states might be the states of maximal entanglement width which can 
be generated via a polynomial number of gates (of fixed squeezing). We leave that comparison for 
future work. 
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to generate states with a greater entanglement-width scaling than for qubit cluster 
states, and therefore that Gaussian cluster states may potentially be rich resources for 
entanglement-based protocols. With the impressive recent progress in the experimental 
generation of Gaussian cluster states [6], developing means to exploit these states for 
information-processing tasks remains a worthy goal, which deserves continued interest. 
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